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1. INTRODUCTION 

We present here the notes of three lectures given by one of us at the XVIII International 
Fall Workshop on Geometry and Physics in Benasque. The course considers some 
aspects of tomographic mappings, a vast and very active area of different applications of 
integral geometry. 

As a motivation, we start with the physical problem of the tomographic imaging of 
a two dimensional body by the measurements of the intensity attenuation of a beam 
of particles passing through it. Then, we will introduce the Radon transform, which is 
the key mathematical tool for reconstructing the tomographic map of both the Wigner 
quasi-distribution of a quantum state and the probability distribution on the phase space 
of a classical particle. The original transform was introduced by J. Radon [22], who 
proved that a differentiable function on the three-dimensional Euclidean space can be 
determined explicitly by means of its integral over the planes. We will consider the 
straightforward generalization of the Radon transform to codimension-one hyperplanes 
in R" and will study its main properties. Moreover, we will prove the original Radon 
inversion formula, by making use of potential theory. 

We will then look at a first generalization of the Radon transform: the M^-transform 
which, roughly speaking, is defined on a larger space. We will show how it can be 
inverted very easily by using the Fourier transform. Then, we will unveil the relation 
between Radon and M^, thus exhibiting an alternative, much simpler inversion formula. 

Then, we will consider a broader framework and look at the more general problem of 
expressing a function on a manifold in terms of its integrals over certain submanifolds. 
This has become an important topic in integral geometry with many applications ranging 
from partial differential equations, group representations, and X-ray technology. The 
focus will be on invariant (or equivariant) transformations under some symmetry groups 
from the space of functions on one geometrical space to the space of functions on 
another geometrical space. We will also exhibit two examples that can be framed in 
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FIGURE 1. Tomographic experimental setup. 



the above setting and work out their respective inversion formulae: the Radon transform 
on codimension-J planes and the John transform on unit spheres in E?. 

We will then show a simple mechanism that allows to generalize the Radon trans- 
form, and its inverse, to more general families of codimensione-one submanifolds. As 
examples, we will look at two families of curves in the plane: the family of all circles 
passing through a fixed point and at the family of all equilateral hyperbolas with one 
fixed asymptote. 

Finally, we will show how the above techniques can be used in the quantum case. 
A straightforward generalization derives from the phase-space description of quantum 
mechanics through Wigner quasi-distribution functions. We will see how the homodyne 
detection scheme can be considered as the Weyl quantization of the Radon transform. 

There is a huge literature on tomography and it is impossible to be complete in 
the references. We apologize in advance for any shortcomings in this respect. There 
are several classical references: we suggest for a group-theoretic approach the book 
by Helgason [9] and for the connection with partial differential equations the book by 
John [12]. Additional general references are Gel'fand and Shilov [8] and Strichartz [23]. 
Finally, for an up-to-date overview of the many applications of the Radon transform we 
refer to the book by Deans [4], which includes also an English translation of the 1917 
paper by Radon, that announced the transform's discovery. 

2. TOMOGRAPHIC IMAGING IN THE PLANE 

Consider a body in the plane M^, and consider a beam of particles (neutrons, electrons. 
X-rays, etc.) emitted by a source. Assume that the initial intensity of the beam is Iq. 
When the particles pass through the body they are absorbed or scattered and the intensity 
of the beam traversing a length As decreases by an amount proportional to the density of 
the body /i, namely 

M=-As^{s), (1) 



so that 

I{s)=IoCxp(^-J\{r)dry (2) 
A detector placed at the exit of the body measures the final intensity I{s), and then from 

= r n{r) dr (3) 
h Jo 

one can record the value of the density integrated on a line. If another ray with a different 
direction is considered, with the same procedure one obtains the value of the integral of 
the density on that line. See Fig 1. 

The mathematical model of the above setup is the following: Given a normalized 
positive and smooth function / = f{q,p), namely f{q,p) > and 

f nq,p)dq&p=\, (4) 

define for all possible lines X in the plane the quantity 

f\X)= [ /(x)dm(x), (5) 

where dm is the Euclidean measure on the line A. In this way, we have defined an 
operator U that maps a smooth function / on the plane into a function on P^, the 
manifold of the lines in M^. Along with the transform f ^ fK it is convenient to consider 
also the dual transform which to a function g on associates the function g" on 

given by 

g\x)^ [ g(A)dju(A), (6) 

JxeX 

where djU is the measure on the compact set (A e P^|;c e A} which is invariant under 
the group of rotations around x and for which the measure of the whole set is 1. 

We ask the following question: If we know the family of integrals ( A ) ) , can we 
reconstruct the total density function /? The answer is affirmative and in the following 
we will see how to obtain this result. 

In order to write the integral (5) in a more convenient form, let us parameterize P^ in 
the following way 

Kx§^p2, (Z,^)^A (7) 

where X G M is the (signed) distance of the line A from the origin, and 
^ = (cos , sin 0) e S is a unit vector perpendicular to A. This means that the Cartesian 
equation of A is 

A = {{q,p) e R^\X-qcosd-psme = 0}. (8) 

By considering a pair (X,y) of Cartesian axes associated to the line A, where X is 
parallel to ^ , and Y is perpendicular to , we obtain the following parametric equation 
for A 

^ = Xcos0-ysin0, /7 = Xsin0 + ycos0. (9) 




See Fig. 2. With a little abuse of notation, we have 

=/B(X,<^) = / dy/(Xcos0-ysin0,Xsin0+ycos0), (10) 

that is called the Radon transform of /. Observe that the (7) is a double covering of P^, 
because the pairs (X,^) and (—X, correspond to the same line A. However, one has 

/«(X,^)=/«(-X,-^) (11) 
and thus depends only on the line A and not on its parameterization. 



3. RADON TRANSFORM 

Let us generalize the above definitions to the case of a w-dimensional space. Let / be a 
function defined on M", integrable on each hyperplane in M" and let be the manifold 
of all hyperplanes in W\ The Radon transform of / is defined as follows 

/»(A) = ^/(x)dm(x), (12) 

where dm is the Euclidean measure on the hyperplane A. Thus we have an operator H, 
the Radon transform, that maps a function / on M" into a function f'^ on P", namely 
f ^ f^. Its dual transform, also called back projection operator, associates to a 

function g on P" the function g'' on M" given by 



g\x)= ^(A)dM(A), (13) 

Jx^X 

where d/j is the unique probability measure on the compact set {A G P^l-^: G A} which 
is invariant under the group of rotations around x. 




Let us consider the following covering of P" 

RxS"-^^r\ (X,<^)h^A, (14) 

where S"^^ is the unit sphere in W\ Thus, the equation of the hyperplane A is 

X = {xeR"\X-^-x = 0}, (15) 

with a ■ b denoting the Euclidean inner product oi a^b & M". See Fig. 3. Observe that 
the pairs (X, {—X, G M x E>"^^ are mapped into the same hyperplane A G P". 
Therefore (14) is a double covering of P". Thus P" has a canonical manifold structure 
with respect to which this covering mapping is differentiable. We thus identify contin- 
uous (differentiable) functions g on P" with continuous (differentiable) functions g on 
Mx§"-1 satisfying ^(X, <^) =^(-X,-<^). 

We will work in the Schwartz space y{R.") of complex- valued rapidly decreasing 
functions on W\ We recall that / G ^(M") if and only if / G C°°(M") and for all 
multiindices a — (ai, . . . , a„) G N" and m > results 



SUp||x|'W/(jc)|<+oo, 

where \x\ denotes the Euclidean x, and 



(16) 



(17) 



In analogy with ^(M") we define ^(M x S"-^) as the space of C°° functions g on 
X S"^^ which for any integers m, a > and any differential operator Z) on S"^^ satisfy 



sup 



IX I 



(18) 



The space ^(P") is then defined as the set of g e ^(M x S"-^) satisfying - ^ ) = 

If / e y{M.") its Radon transform (12) can be written in a simple way, by using the 
Dirac 5 distribution, as 

fix)d{X-^-x)dx^{d{X-^.x))f, (19) 

Jw 

In the same way if g e y{¥") one has 

g\x)= f g{X,^)d{X-^-x)dXd^. (20) 

Using (19) and (20) it is easy to check the duality relation between the Radon transform 
tl and its dual b, namely for every / e ^(R") and g e ^(P") one gets 

[ dXd^/»(X,Og(X,0= / d^/W/W- (21) 

Indeed, 

/ 6Xd^f\x,^)g{X,^) = [ dXdU dxf{x)5{X-^.x)g{X,^) 

dx/(x)g^(x). (22) 



Now, consider a unit vector e ^ We define tomogram of f along t, the 
function 

/|(X):=/«(Z,0. (23) 

Physically, this means that we fix a direction ^ and we look at the variation of the 
intensity of all the rays that interact with the body in that direction. We also define the 
family (/|)^£§«-i of tomograms associated to /. 

The following result shows how the Radon transform maps a joint probability den- 
sity function f of n random variables into a family of probability density functions 

(/|)^eS«-i °f ^ random variable. 

Lemma 1. Let / G ^(M") be a nonnegative and normalized function, namely 

f>0 and / /(x)dx=l. (24) 
Jw 

Then, for all e S"~^ it results that 

fl>0 and [ fUx)dX=l. (25) 

^ Jr ^ 

Proof Let us fix ^ e §«-^ Observe that for all Z e M 

/ f{x)d{X-^-x)dx>0. (26) 



Moreover, 

[ fl{X)dX= [ dX [ f{x)d{X-^-x)dx^ [ f{x)dx=l. (27) 

□ 

The Radon transform has a pecuUar behavior under translation. Given a displacement 
a e R", consider the translation operator defined by 

{Taf){x)=f{x-a). (28) 

One has the following property 

Lemma 2. Let f e y{W) and a e R". Then for any ^ e §""1 

M)l=^.afl (29) 

Proof. By setting y = x — a we get 

{TafiiX) = / f{x-a)d{X-^.x)dx^ [ f{y)d{X-^.y-^.a)dy 

□ 

The following lemma, which is the infinitesimal version of Lemma 2, shows how 
the Radon transform and its dual intertwine the Laplacian on R", A = d^. with the 
Laplacian on R, 

Lemma 3. Let f e y{R") and g e ^(P"). Then 

(A/)« = ^ and (^^J ^^i. (30) 
Proof. By Lemma 2 we get 

(T«/)»(X,0=/"(^-^«,0- (31) 

Therefore, by considering an orthonormal basis (ei, . . . , e„) of R", for all X e R" and 
^ e§""^ one has that 



and thus 

dx} 



Therefore, since G §" \ by summing over i we have the first equaUty in (30). 
As for the second relation, observe that 

g\x)=[ g{X,^)5{X-^.x)dXd^= [ e(^x,^)d^, (34) 

whence 
Therefore, 

VW = £ ^ W - |^(« d« = W. (36) 

□ 

The Radon transform is closely connected with the Fourier transform. Let us recall 
some basic facts about the Fourier transform. For any function g e ^(M™) with ineN^, 
we define its Fourier transform as 

g{k)= I gWe-'*-"ic, (37) 

with k e W^. It is not difficult to prove that the Fourier transform is a bijection of S^{W^) 
onto itself, and its inverse reads 

g{x)= I g{k)e'-^j^. (38) 

In short, for any g G ,5^{W"), one has that g = g = g. Using standard duality arguments 
the Fourier transform can be extended to the space of temperated distributions y'iW^) 
(and thus also to its subspace of square summable functions L^(]R"), where it can also 
be proved to be unitary). In particular, we recall that 5 = 1 and thus 

5(x)= I e'^-^-^. (39) 

We will show that the n-dimensional Fourier transform of a function / on M" is the 1- 
dimensional Fourier transform of its tomogram /| for all ^ G S"~^ This result is known 
as the Fourier slice theorem. 

Lemma 4. Consider f e ^{W), t G K an^i ^ G S''-^. Then 

fiz^)^fl{z) (40) 



Proof. The thesis follows immediately from the definitions. Indeed 

fUx) = [ dX c-''^ fUx) ^ [ dX e^''^ [ dxf{x)d{X-^-x) 

^ JR ^ JR JR" 

= [ dxfix) [ dXe-'^^5{X-^-x)= [ dx/(x)e-'^^-^ = /(T<^). (41) 

jR" JR JR" 

□ 

Now we will address the problem of the range of the Radon transform with domain 

o$^(]R"). The following lemma is crucial. 

Lemma 5. /// e y{R"), then for any integer A; e N and any ^ = (^i, . . . , e §"~^ 
the integral of its tomogram, 

If{^)= [ fliX)X'^dX, (42) 

JR 

is a homogeneous polynomial of degree k in the Cartesian components ^i, . . . , 
Proof. This is immediate from the relation 

[ f\X,^)X'^dX^ [ dXX' [ dxfix)d{X-^-x)^ [ f{x)i^-xfdx. (43) 
JR Jr JR" JR" 

□ 

In accordance with this lemma we define the space as follows: g e J^(P") 

if and only if g G J^(P") and for any integer A: e N the integral 

7,(0= / 8iX,^)x'dX (44) 

JR 

is a homogeneous polynomial of degree km^i,...,^n- 

Theorem 1. The Radon transform is a linear one-to-one mapping of y(M.") onto 
^h(P«). 

The proof uses the well-known fact that the Fourier transform maps onto itself 

and Lemma 5. For more details see [9]. In the next section we will give a constructive 
proof, by explicitly exhibiting the inverse mapping. 

4. INVERSION FORMULA 

Now we want to obtain the inversion formula, namely we want to prove that one can 
recover a function / on R" if all its tomograms (/|)^£§n-i are known. In order to get 
this result we need three preliminary lemmata. 

Lemma 6. Let f e and V(x) = 1 /\x\, xeM.",xj^O. Then 

{fi)\x)=anf*V, (45) 



potential {P)\x) 




charge distribution f{y) 

FIGURE 4. {p)^{x) is the potential at x generated by the charge distribution /. 

where a„ depends only on the dimension n, and * denotes the convolution product, 

(/*^)W= / ny)gix-y)dy. (46) 

Proof. Observe that 

(/¥W = / d^M-x,^)=f dyf{y)5i^-x-^-y) 

= [ dU dyfiy) [ ^e-^ (^--^). (47) 
j§"-i Jr" jr 2k 

Now, with the substitution x — 3; =\x — y\o),G)E and t = s\x — y\,ihQ integral (47) 
becomes 

/ d^f dyny)[^y^e^'^"^=a,J dy fiy)Vix-y), (48) 
J§"-i JR" jr2h \x — y\ JR" 

where 

an= [ ^[ d^e^^^-. (49) 

By symmetry, it is apparent that a„ is independent of CO. □ 

The physical meaning of Lemma 6 is the following: if / is a charge distribution, 
then the potential at the point x generated by that charge is exactly {f'^)^{x), see Fig. 4. 
Notice, however, that the potential of a point charge scales always as the inverse distance 
independently of the dimension n, and thus it is Coulomb only for n = 3. The only 
dependence on n is in the strength of the elementary charge a„. This fact will be crucial in 
the following: we will see that the associated Poisson equation involves an n-dependent 
(fractional) power of the Laplacian. 

Lemma 7. Let V{x) = l/\x\, x & M", x^O . Then its Fourier transform is 



(50) 



where bn depends only on n 
Proof. One gets 

Vik)^ [ d;cAe-i^-=/ dx^c-'\'\''-^= [ -^le-i«-^ = -^, (51) 

JR" \x\ jRn \x\ Jw i \y\ \k\"' ^ 

where we have set k— \k\co with co e S"~\ y — \k\x, and 

bn = / dx^— . (52) 

JR« \X\ 

Moreover, bn is independent of a by symmetry. □ 

Lemma 8. Let an and bn as in Lemmata 6 and 7. Then 

anbn = 2{2TtY-\ (53) 

Proof. Since the Fourier transform is a bijection, we use the trick of going back and 
forth. From Lemma 7 we have 

V(x) = ^ = ^{x)= [ Vik)e'-^-^ = [ -^ei*--^. (54) 

By setting = |j:|(0 and p = \x\k = t^, with CO, ^ G S"~^ and t e M+, Eq. (54) becomes 

]_ ^ bn j dk ^ifc.^i^i ^ bn j dp ^ip.^ 



(2;r)"|x| ,/e+ Js"-! 



d£ e^^''" = " " ^ (55) 

2(27r)"-i Je 2;r i§«-i ^ 2(27r)«-i |jc| ' 

where in the last equality we used Eq. (49) of Lemma 6. Therefore, 

, = 1 (56) 
2(2;r)«-i 

and the thesis follows. □ 

Now we are ready to prove the inversion formula for the Radon transform. 
Theorem 2. Let f e ^(R"). Then 

/W = ^(^(-A)'^(/VW (57) 
where {—A)", with a >0, is a pseudodifferential operator whose action is 

((-A)«/)W=/ |^P"/»e''-^7:^. (58) 



Proof. First of all observe that if n = 1 statement (62) is trivial, so assume that n> \. 
Let p > 0, we get by Lemma 6 



{-Arv{x)^[ j^Jk\'PV{k)c^-^ = bJ 7^|^P^-"+^ei^ (59) 
so that, if 2/7 = n — 1, we have 

(-A)'^V(;c) = MW. (60) 
Therefore, again by Lemma 6, we have 

(_A)¥(/«)b(^) = a„ (/*(-A)'2^y) {x)^anbn{f*d){x)^a„bnf{x). (61) 

From (53) we get the inversion formula (62). □ 

n— 1 

Equation (57) says that, modulo the final action of (— A)t^, the function / can be 
recovered from its Radon transform f'^ by the application of the dual mapping b: first 
one integrates over the set of points in a hyperplane and then one integrates over the set 
of hyperplanes passing through a given point. Explicitly we get 

which has the following remarkable interpretation. Note that if one fixes a direction 

^ G then the function f{^-x,^) is a plane wave, i.e. a function constant on each 
plane perpendicular to ^ . Therefore, Eq. (62) gives a representation of / in terms of 
a continuous superposition of plane waves. A well-known analogous decomposition is 
given by Fourier transform. When n = 3, one recovers the inversion formula originally 
found by Radon [22] 

/W = -^A/^^/«(^x,^)d^ (63) 
4.1. Localization 

Now we consider the following question on the localization properties of the Radon 

transform: given a function / can we reconstruct it in a compact region Q., using only 
the tomograms of / on hyperplanes A that pass through D.1 By looking at the inversion 
formula (62) one would naively think that the answer be positive. By following a nice 
argument given by John [12], we will show that in fact this cannot be true in two (and, 
in general, even) dimensions. The reason is the well-known fact that waves cannot be 
localized in two dimensions: a pebble dropped in a pond produces ripples that propagate 
along expanding disks. On the other hand, in one dimension waves remain localized: a 
local perturbation of the pressure in an organ pipe propagates by translating its support. 
This is a fortunate accident that allows us to enjoy a Bach's fugue without hearing an 
endless annoying echo of each tune. 



FIGURE 5. A function / supported in a ball and its tomogram in the direction ^ . 



Consider the following Cauchy problem for the wave equation in M?: 

d^f-Af = 0, /(0,x) = g{x), drmx) = 0, (64) 

and assume that the initial condition g has its support in the ball 5(0, £) = { |x| < e} with 
center and radius £ > 0. Let us fix a time t > e. The solution to this equation is given 
by the celebrated Poisson's formula [5] 

fM = ^[ ^MlSkH^d,, (65) 

where V denotes the gradient in M?. Therefore the support of /(?,■) is a subset of 
B{Oj + e). On the other hand, if we consider the Radon transform of (64) we get, for 
every ^ G ^ , a Cauchy problem for the wave equation in M 

dlfl-dtl=0. /«(0,X)=4(X), aj«(0,X)=0. (66) 
Observe that has its support in (— e, e). The solution of Eq. (66) is 

/«(f,X) = i +gliX + t)) , (67) 

thus the function /|(?, ■) has its support in (—? — £, —? + e) U (? — £,? + £). So we see that 

if one wants to reconstruct the function /(?, ■) in the set 5(0, f — e) one cannot consider 
only the tomograms on the lines that intersect 5(0,? — e), because they are all equal to 
zero. See Fig. 5. Indeed, from f^{t,X,^) = for \X\ <t — £ we have that 

(/«)^(?,x)= //tt(r,^x,^)d^=0, for \x\<t~£. (68) 



On the other hand, the inversion formula (62) forn = 2 yields 

nt,x) = j^{-Ay/\fi)\t,x)j^O, for \x\<t-e. (69) 

The seeming contradiction derives from the fact that the pseudodifferential operator 
(— A)^/^ is nonlocal. In fact, by (58) it is an integral operator that depends also on the 
values of {f'^)\t,x) for \x\ > t — e. Of course, this is a peculiar behavior of the even 
dimension. For odd n in the inversion formula (62) we have an integer power of the 
Laplacian that is a bona fide differential operator and thus is local. 



5. M^-TRANSFORM 



In this section we introduce another covering of the homogeneous space P" that uses an 
additional parameter, which is absent in the classical tomographic framework of Radon. 
This choice of parametrization of P" gives an easier inversion formula, although this 
approach is completely equivalent to that considered in the previous section. 
We define the following covering 

MxM'J^P^ (X,il)^X, (70) 

where the Cartesian equation for X is 

A = {xeR"|X-;U-x = 0}. (71) 

Of course the additional parameter with respect to Radon is |ju|. Now, if / e we 
can define the -transform of / as follows 

f^\x,ii)^[ f{x)d{X-ii-x)dx^{d{X-ii-x))f. (72) 
Jr" 

This tomographic mapping was introduced by Man'ko and Marmo, who gave seminal 
contribution towards its significance [16, 17, 18]. Also for this transform we can define 

the tomograms, precisely for every ;U G M" we define f^f {X) :— It easy 

to check that also this transform maps a probability density on R" into a family of 
probability density functions M, namely if / e y{M."') is a positive and normalized 

function, then for all jU G M^, fjf is positive and normalized. 

Let us now consider the relations between Radon and transform. Obviously, from 

f'^ {X,^.) one can immediately recover f'^{X,^)hy setting jU = ^ G §"~^ : 

/tt(Z,^) = /^'(X,.^). (73) 

However, notice that, although /* is the restriction of f^^ to the unit sphere S"^ ^ , there is 
actually a bijection between the two transforms. Therefore, they carry exactly the same 
information. Indeed, since the Dirac distribution is positive homogeneous of degree — 1, 
i.e. 5 (ax) = \a\~^5{x),fov every a ^ 0, one gets from Eq. (72) 

(74) 



for /i 7^ 0. In words, the tomogram f^^(X,[i) at a generic point /i G M" is completely 
determined by the tomogram at /i/|ju| G S"~^ But the latter is nothing but the Radon 
transform, by Eq. (73). Therefore, we get the bijection 

Ax,^)=f^\x,^), f^\x,^) = . (75) 

IMI \ IMI IMI / 

Finally, let us look at an inversion formula as in the case of the Radon transform. 
Theorem 3. Let f G ^(R"). Then 

Proof. Observe that 

/ Axf{x)5{X-^.x)= f Axf{x) f ^ei^(^-^-) 

An+i (2;r)«+i -^Uy |t|« L+i (271)"+^ ^ 8W),K//) 

where g{T,y) = {2nYf{y/x)/\x\''. Therefore, /^^(X,ju) = g(X, -ju), and so 

g(T,x) = P(T,-x). (78) 

Finally, 

1,-x) _ f^{-\,x) _ f dXdii 



Summarizing, the inversion formulae for the transforms (62) and (76) read 



(79) 

□ 



/W = / /^'(X,M)e'(^-^-^)^, (81) 

respectively. While formula (81), which is nothing but a Fourier transform, is quite easy 
to handle, the inversion formula (80) is in general very hard to tackle, especially for even 
n, due to the presence of a fractional Laplacian. For an explicit example, see [6]. 



6. GROUP-THEORETIC FRAMEWORK 



The general problem posed by the Radon transform is that of determining a function on a 
manifold by means of its integrals over a certain family of submanifolds. Here, following 
Helgason [9] we will construct a quite general framework for the above problem. 

Let us start from the simple observation that the Euclidean group E{n) acts transitively 
in a natural way both on the affine space R" and on the hyperplane space P". The 
Euclidean group is the group of isometrics associated with the Euclidean metric and 
is the semidirect product of the orthogonal group 0{n) of rotations, extended by the 
abelian group of translations M", namely E{n) — 0{n) kW. 

Fix an origin in R" (a reference point) and an origin Aq in P" (a reference hyper- 
plane). It will be convenient to choose an hyperplane Aq passing through the origin 0. 
Transitivity means that any element of R" (P") can be obtained by an Euclidean motion 
of the origin (Aq), namely 

\/x eR" : x = g-0, for some g 

VA e P" : A = 7- Ao, for some reE{n), (82) 

where the dot • denotes the action of E{n). The action is not free, i.e. one can have 
g ■x = h-x for some x € R" and g,h E E{n) with g ^ h. Therefore, let us consider the 
isotropy subgroup of the origin (also called the stabilizer subgroup) as the set of all 
elements in E{n) that fix 0: 

E{n)o = {ge E{n) | g • = 0} = 0{n). (83) 

It is given by the rotation group (around 0). On the other hand, the isotropy subgroup of 
the reference hyperplane Aq passing through the origin is 

E{n)^ = {ge E{n) | g ■ Ao = Ao} = Z2 x - 1), (84) 

where E{n — 1) is the Euclidean group of the (n — 1) -dimensional affine subspace Aq 
(translations and rotations in Aq), while Z2 ~ 0(1) is the group whose elements are 
the identity and the reflection in Aq. See Fig. 6. Now, given an element x e R", by 
transitivity we get x = g - 0, for some g E E{n). Therefore we can identify x with 
the left coset gE{n)Q = gO{n). Analogously, for all A G P" we get the identification 
A ~ yE{n)x^ — 7(Z2 x £(« — 1)), for some y EE{n). Therefore, we get that 

R" ~ E{n)/E{n)Q = E{n)/0{n), 

P" ~ E{n)/E{n))^=E{n)/{Z2xE{n-\)), (85) 

are homogeneous spaces of E{n). 

We want to describe the incidence relation in this framework. Namely, given an 
hyperplane A G P" and a point x G M", when x lies in A, that is x G A? It is not difficult 
to see that this is true if and only if their corresponding cosets, considered as subsets of 
E{n), have a point in common. In equations, given x ~ gO{n) and A ~ 7(Z2 xE{n—\)), 
one has 

xeX ^ gO(n)n7(Z2x£(n-l))^0. (86) 




Indeed, one gets 

g-Oey^ ^ g-0 = Yh-0, for some /? e Z2 X ^(n- 1) 
gk = yh, for some k E 0{n) 
^ gO{n)r}y{Z2xE{n-\)) = {gk}. (87) 

In the Cartesian product M" x P" we can consider the following set, called the inci- 
dence relation, 

/={(x,A) gM"x P"|xe A}. (88) 

Consider now the restrictions to / of the natural projections p : / C M" x P" — )• M" and 
;r : / C X P" — > P". Obviously, since every point of M" lies in some hyperplane of P" 
and every hyperplane passes through some point, we get p{I) = W and = P", so 
that p and n are onto. Now, for any x G and A G P" consider the sets 

x^ = {XeF"\{x,X)el}, X^ = {xeR''\ix,X)el}. (89) 

They represent the set of all hyperplanes passing through x and the set of all points lying 
in A, respectively. The crucial fact is that the mappings 

x^x^ = 7rop-\x), X^X^=po7j:^\X), (90) 

are in fact injective. Indeed, we obviously have 

x= f]X, X= \Jx. (91) 

Therefore, any point x G M" can be uniquely identified with a set of hyperplanes x' C P" 
(those incident to x), and conversely any hyperplane X G P" can be uniquely identified 
with a set of points A'* C M" (those incident to A). We say that M" and P" are homoge- 
neous spaces in duality. In this setting we can rewrite the Radon transform and its dual 



as follows 

fK^) = /^/W dmW, q>\x) = (p(A) dM(A). (92) 
Notice that from Eq. (87) we get for any pair (x, A) e / 

(x, A) ~ {gO{n),gk{Z2 X £(n - 1))) = {glO{n),gkl{Z2 x £(n - 1))) , (93) 

for some k e 0{n) and for all / G L, where 

L = E{n)Qr\E{n))^ = 0{n) n (Z2 x £(n - 1)) = Z2 x 0(n - 1) (94) 

is the subgroup that fixes both and Aq. Therefore, we get that also the incidence relation 
7 is a homogeneous space, namely 

I^E{n)/L = E{n)/{Z2xO{n-\)). (95) 

The Radon transform on the homogeneous spaces in duality R" and P" can be conve- 
niently described by the following double fibration: 

7~£(n)/(Z2xC>(n-l)) 



M«~£(n)/0(n) P"~£(n)/(Z2x£(n-l)) 

(96) 

Now let us generalize the above construction. Let G be a locally compact group, let K 
and H be two closed subgroups of G, and \q\.L = KHH. Assume that KH C G is closed. 
Let us consider the following two left coset spaces 

X = G/K and A = G/H, (97) 
and define the incidence relation 

7 = {(jc,A) eXxA|jcnA 7^0}. (98) 
Define also for all x e Z and A e A 

x^ = :;rop-i(x) = {A e A| (x,A) e/}, = pon'^X) = {x eX\ {x,X) e 1} , (99) 

where p and 7t are the restrictions to I of the natural projections of Z x A. It is not 
difficult to prove thatx^ andA^ are closed subsets, and that the sets 

Kh^ {keK\kHUk^^H cHK} , Hk = {h e H\ hKUh^^K C KH} (100) 

are subgroups of K and H, respectively. Moreover, Helgason [9] proved the following 
Theorem 4. The following propositions are equivalent: 
1. 

L = Kh = Hk (101) 



2. The mappings 



(102) 



are injective. 

Therefore, under the trans vers ality assumption (101), any element x of X can be 
uniquely identified with the subset }^ of A, and any element A G A can be uniquely 
identified with the subset C X. Thus, X and A are homogeneous spaces in duality. 
Exactly as in the case of hyperplanes in R", we have the identification I = G/L, with 
L = KDH, that is easily proved via the bijection G/L ^ I : gL^ {gK,gH). As a 
consequence, also in the general case we have the following double fibration, which 
is the generalization of (96): 



I = G/{KnH) 




(103) 



X = G/K A = G/H 

Let xq and Xq be the origins of X and A, respectively. Let us assume that the sets Xq = 

K/L and Xq=H/L have positive measures invariant under K and H, respectively. Under 

the transversality assumption (101), by using the action of G one can uniquely induce 
a ^-invariant (//-invariant) positive measure /i (m) on each x^ (A^*). Therefore, one can 
define a Radon transform and its dual as in Eq. (92) for the double fibration (103). 

A final remark is in order. The above construction in terms of homogeneous spaces in 
duality is quite suggestive and represents a good framework for encompassing various 
examples of tomographic mappings. However, as noted by Helgason himself [9], it is 
doubtful that it could capture the large variety of results in this field. In particular, there 
remain several open problems that should be addressed case by case. In particular, one 
should find suitable function spaces on X and A where the Radon transform / i-> /l* is 
invertible, and - an even more difficult task - find an inversion formula. 

7. EXAMPLES 

In this section we consider two examples of tomographics mappings a functions/ on R", 
in terms of its integral over certain submanifolds of M" and prove the related inversion 
formulae. In the first example we consider the planes of codimension d, < d < n, in 
M", while in the second one we consider the unit spheres in M?. 



7.1. (n-d)-planes in M" 

Let us consider the space of all the affine subspaces A of M" with dimension n — d, 
< d < n.la this example the homogeneous spaces in duality are R" and and the 



double fibration in Eq. (103) reads 



I = E{n)/{0{d) X 0{n-d)) 




M" = E{n) join) = E{n) /{0{d) xE{n-d)) 

(104) 

Indeed, G = E{n) and K = 0{n), the isotropy group of the origin 0, hence X = R" = 
E{n) join) . Moreover, the isotropy group of a reference subspace Ao G passing trough 
theoriginis// = C>(J) x£(n- J), and thus A = P^ = x if)). Finally, 

L^KHH ^0{d)xO{n-d), so ±SLtI ^ E{n)/{0{d) xO{n-d)). 

The Radon transform has the same form (12) of the codimension-one case, where 
dm{x) is the Euclidean measure on the codimension-if hyperplane A. First of all notice 
that, since a codimension-one hyperplane can be viewed as a disjoint union of parallel 
codimension-J planes, it is obvious that the transform / /* is injective. In other 
words, the information recorded is redundant and we can recover the function / by 
integrating away d—l degrees of freedom and then inverting the usual Radon transform 
on hyperplanes. We will follow here a more direct, different approach. 

Observe that A can be obtained by the intersection of d orthogonal hyperplanes in M", 
namely A can be described as the space of the solutions of the following linear system 

X = {xeR"\Xj = ^j-x, 1,...,J}. (105) 

where Xj G IR and G E>"^\ ■ = Sj.k, with j,k = 1, ... .J. Of course, this repre- 
sentation of A is not unique, in fact we can obtain A as the solution of another linear 
system 

A = {x e I F^- = rij -x, 7 = 1, . . . ,d}. (106) 
where for all / = 1, . . . , J 

d d 
li=I,ciij^j, Yi=Y,aijXj, (107) 

;=i 7=1 

and A = (o/y) G GL(J,R). Let us compute 

d 

k,l=l 

Note that, since rj, G S"^\ one must have that (AA^),,- = 1 for all z = 1, . . . 

Let / G ^(R") and define the tomographic transform of / on A G P^ as follows (with 
the usual abuse of notation) 

/tt(A) = /tt(F,T])= / fix) n - r]j-x) Vdet(T]^ • m) dx, (109) 

Jr" jj[ 



where F = (Fi , . . . , F^) G and r] = (r]i , . . . , 7]^) G {S"-^Y. The definition of the Radon 
transform in (109) does not depend on the parametrization of A. Indeed, by writing 
the parameterization (F, 77 ) in terms of a parameterization with orthogonal hyperplanes 
(X, (^) through (107), we get 



dt 

d 



^ [ ^iLjij-LiajiiXi-^rx) idetAl = / e^^'fe''^''^)^^'-^'-^) IdetAl 

= / ei^'^'(M-) = n5(Z,-^;.x). (110) 

where we set s = A^t. Therefore, 
f^{Y,rj)=f{X,^)= [ f{x)fl8{Xj-^j-x)dx= [ f{x)dm{x)=f^{X), (111) 

as given in (12). 

Let us consider also the dual of the Radon transform. If g G ^(P^), we define 

8\^)-l ^ ^^AY,n){[d{Yj-nj.x)^£^= (112) 
J{RxS"-^Y AA ydet(T];t-77/) 

Theorem 5. Let f e ^(M"). Then 

/(.) = c„,,(-A)^(/«)^(.), ^here = ^ jixl) J/^L^ , fii^) 
w«Y/i r denoting the Euler gamma function. 

Proof. First of all observe that if n = 1 then necessarily d = \, and the statement (113) 
is trivial. So assume that n> I. Let < a < n, then if y(x) — l/\x\, with x G R", one 
can easily prove that [14] 



V«-« = Z7„,«y«, where = 2«;r«/2-^Ali^ 

^ V 2 j 



(114) 



Let p>0 and let us compute the following quantity 

dk 



i-Arv'ix) = / T^r-^ |^|2^^v5(^)e*- (115) 
SO if 2/7 = n — J, we obtain 

{-Ay-^V^{x) = b„,„.dS{x). (116) 



Now let us compose the Radon transform with its dual. As in Lemma 6 we obtain 
(/.fW=.^(/.V^)W, ., = /J|_,,|e»^«. = ^ = ^. (1,7, 

and again a„ is independent of ft) e S"~^. Therefore, we have 

{-^)"^{fi> = ai{-A)"-^{f.V')=aif.{{-A)"-^V') 

= aiK,n-df*^^aiK,n-df (118) 
and so we have Eq. (113) with c^^ = l/{a^bn^n-d)'> this concludes the proof. □ 

7.2. Spherical means in 

Let us consider the case, studied by John [12], of a tomographic transform that maps a 
function / on M" into a function on the spheres in M" with a fixed radius r. For simplicity 
we consider the case n = 3 and r = 1. 

The homogeneous spaces in duality, and = set of all unit spheres in 
(parametrized by their centers), are described by the following double fibration: 



I = E{3)/0{2) 




(119) 



R^ = £:(3)/0(3) m^=E{3)/d{3) 

Indeed, G = E{3) and K = 0(3), the isotropy group of the origin 0, hence X = = 
£■(3)70(3). Moreover, the isotropy group of a reference unit sphere Ao G passing 
trough the origin is H = (5(3) (rotation around the center of Aq), and thus A = R^ = 
E{3) /d{3). Finally, L = KnH = 0(3) n 0(3) = 0(2) (rotation around the line joining 
and the center of Ao), so that / = £(3) /0(2). 

We define the tomogram of / as its average over the unit sphere with center x 

f^i^)=lL[ji^+^)d^^ (120) 

1^ I J§2 

where | = 4;r is the surface of S^. We look for an inversion formula, namely we want 
to reconstruct the function / using its tomograms on all possible unit spheres in R^. This 
problem is equivalent to the following integral equation for the function /: 

±-J^J{x+^)d^=g{x). (121) 



In Ref. [12] John proved that the solution of the integral equation (121) is unique, under 
mild regularity assumptions on / and g. Therefore, also this tomographic mapping is one 




FIGURE 7. Inversion formula (128). Dotted lines are spheres with an odd radius, while dashed lines 
are spheres with an even radius. 

to one on y{M?). We will now explicitly construct its inverse. To this purpose, consider 
the following Cauchy problem for the wave equation 

d^v-Av = 0, v(jc,0)=0, dtv{x,0)=g{x). (122) 

The solution of this problem is 

vix,t)=tl^s\x,t), iis)(x,t) = r^ [ g{x + t^)d^. (123) 

We construct the function u{xj) as follows 

u{x,t) = - [dtv{x,2k+l+t)-dtv{x,2k+l-t)] (124) 

and define f{x) := dtu{x,0). Observe that u{x,0) = 0, thus u solves the Cauchy problem 
d^u-Au = 0, m(x,0)=0, dtu{x,0) = f{x), (125) 
whence u{x,t) = tl^f\x,t). In particular, we have 

-\-oo 

u{x, 1) = 1) = - £ [dtv{x,2k + 2) - dtv{x,2k)] = dtv{x,0) = g{x) (126) 

and so we have proved that the function / solves Eq. (121). Moreover, 

-|-oo 

/(jc) = a,M(;c,0) = -2£^2^(x,2A:+l) 

k=Q 

-\-oo -\-oo 

= -2£ Av(x,2yt+1) = -2AY^{2k+l)I^^\x,2k+l) 

= -2^1-^/ dT]/(x + (2fc+l)^+7]). (127) 




Therefore, we finally obtain 

-\-oo r\ 1 I 1 p 

fix) = -2A£ — ^ d^f^{x+{2k+l)^), (128) 

k=o PI •^s^ 

which is a sum of the spherical means over all the spheres centered at x and with odd 
radius. See Fig. 7. 



8. TOMOGRAMS ON HYPERSURFACES 



A simple mechanism that allows nonlinear generalizations of the Radon transform is 
the combination of the standard transform with a diffeomorphism of the underlying 
M" space [2]. Let us consider a function f{q) on the n-dimensional space q eR". The 
problem is to reconstruct / from its integrals over an n-parameter family of submanifolds 
of codimension one. 

We can construct such a family by diffeomorphic deformations of the hyperplanes (in 
the X eM" space) 

X-ii-x = 0, (129) 
with X G M and ji G M". Let us consider a diffeomorphism of M" 

qER" ^x=(p{q) ER". (130) 

The hyperplanes (129) are deformed by (p into a family of submanifolds (in the q space) 

X-ii-(p{q)=0. (131) 

The case n = 2 is displayed in Fig. 8. 

Given a probability density f{x) on the x space, the -transform can be rewritten as 



where 



f{x) d{X - iJ.-x)dx 
J{q) 



fi(piq))d{X-II-(p{q))J{q)dq, (132) 



(133) 



dxi 




d(pi{q) 


dqj 




dqj 



is the Jacobian of the transformation. 

Observe now that f{x) dx — f{(p{q))J{q) dq, whence 

M = K<(>{q))J{q) (i34) 

is a probability density. Therefore the tomograms along the codimension-1 submani- 
folds (131) are given by 

f'P{X,ll)^{5{X-ll-(p{q)))f^ [ f{q)5{X-ll-(p{q))dq, (135) 

withX eRandjU eM'J. 

The inverse transform follows by (76): 

f{q)^f{cp{q))m = /^^^/n^,M)-/(?)e^(^-^-^(^)) (136) 



with a modified kernel 



dqj 



(137) 



Therefore, a probability density distribution on M" 

/(?)>0, [ f{q)dq=l, (138) 
produces tomograms f^{X,ii) that are probability densities 

f'P{X,ii)>0, //'P(X,ju)dX = l, WiieK. (139) 

JR 

We will now consider two applications of these generalizations of the Radon transform. 
For additional examples see [1, 2]. 

8.1. Circles in the plane 

In the punctured (x^y) plane without the origin (0, 0), the conformal inversion 

{x,y) = cp{q,p) = f -2^) , (140) 

maps the family of lines X — jux — Vy = into a family of circles 

X{q^ + p^)-^q-vp^0, (141) 

centered at C = {jl/lX, v/2X) and passing through the origin (see Fig. 9). When X = 
they degenerate into lines through the origin. The Jacobian reads 

J{q,p) = {q^+pY\ (142) 




whence the transformation is a diffeomorphism of the punctured plane onto itself. 
The singularity of the transformation at the origin (0, 0) is irrelevant for tomographic 
integral transforms because it only affects a zero measure set. 
Equations (135)-(136) become 

r{X,li,v) = l5(x-^^-^^)) , (143) 
f{q,p) = / jnX,^) dXd^dv. (144) 



8.2. Hyperbolas in the plane 

In the {x,y) plane the family of lines X — jUx — = is mapped into a family of 
hyperbolas 

X-^-Vp = 0, (145) 

q 

with asymptotes q = and p = X/y , by the transformation 

{x,y) = (p{q,p) = {q-\p). (146) 

For II > the hyperbolas are in the second and fourth quadrants, while for /i < they 
are in the first and third quadrants (see Fig. 9). When /i = or v = they degenerate 
into horizontal or vertical lines, respectively. The Jacobian reads 

J{q,p)=q-^ (147) 

whence the transformation is a diffeomorphism in the cut plane without the axis (0, j). 



Equations (135)-(136) become 



f<P(X,n,v)^(^d(x-^-vpj^^, (148) 

9. QUANTUM TOMOGRAPHY 

We now consider a possible generalization of the above classical techniques to the 
quantum case. The most direct path is that of using the phase-space description of 
quantum mechanics through Wigner quasi-distribution functions [25, 20, 10]. To this 
purpose, let us consider the phase space r*R'' ~ R^" of a classical system with n 
degrees of freedom. Let us assume that the (statistical) state of the system is given by a 
probability density function f{q,p) in The expectation value of an observable, 

which is represented by a function o{q,p) on the phase space, in the state / is given by 

(a)/ - / a{q,p)f{q,p)dqdp. (150) 
The -transform (72) of the density /, and its inverse (76) read 

/^'(X,AI,V) = {d{X-ll-q-V-p))f= [ f{q,p)d{X-II-q-V-p)dqdp, 
fM = 1 /^^^,M,v)eK^-^-^-^)'^. (151) 

with Z e R and (/i, v) e R^". Note that the tomogram /^J^^y) is nothing but the prob- 
ability density function of the observable X{q,p) = jj, ■ q+v ■ p. A very special set 
of observables are the canonical coordinates: position q = (<?i, . . . ,qn) and momentum 
p = (pi, . . . ,Pn)- One obtains their density directly from the tomograms (151) along the 
Cartesian axes. At jU, = 5ij and v = one has X{q,p) — qj and thus 

/^'(X,(5y),0) = J^J{q,p)d{X-qj) dqdp 

= J^^^ Jiiqi, ■ ■ ■ ,qj-l,X,qj+i, . . . ,qn),p) dqi . . .dqj-idqj+i . . .dq„dp, (152) 

that is the marginal density of the position component qj. Analogously, at jU = and 
V, = 5ij one has X{q,p) — pj and 

/^'(Z,0,(5y))= / f{q,p)d{X-pj)dqdp (153) 



is the marginal density of pj. 



In quantum mechanics (more precisely, in its Schrodinger representation) the Carte- 
sian canonical coordinates q and p are promoted to operators, q and p, that act in the 
Hilbert space of square summable functions L^(]R") endowed with the scalar product 

(V^,0)= / v^(x)0(x)dx, (154) 

with V'^, ^ G L^(]R") and the overline denoting complex conjugation. Precisely, q is the 
multiplication operator {lq^f){x) — x^f{x), and the p — —ihV is a differential operator, 
both with maximal domains. Here, h is the Planck constant. Pure quantum states of 
the system are associated to normalized elements e L^(R"), called wave functions, 
II V^lli ~ {W^W) = 1- The expectation value of position (momentum) in the wave function 
i/A is given by the sandwich (i//, ^i//) ((i//,pi//)). It easily follows that x i— )> |i//(jc)p is 
the marginal probability density of position (correctly normalized, since ||v^||2 = !)• 
Moreover, by Fourier transform it is also not difficult to show that the marginal of 
momentum is given by p i— )> (IkH)" (prove it!). 

Now we want to define the quantization of any classical observable o (also called 
symbol), that for simplicity we assume in Observe that 

aiq,p) = I d{n,y)S^-^^yp) (155) 

SO, by following Weyl, we can construct an operator out of o by substituting q with the 
position operator q and p with the momentum operator p. Thus, we are led to consider 
the unitary operators 

(y, T]) ^ (e^in-9+yp)yf^ = Q^nix+hy/2)^(^^_^j^y^^ ^^5^^ 

(the Schrodinger representation on L^(]R") of the Heisenberg group [7]), with {y,r}) G 
M^" and \j/eL^ (R") . In conclusion, for any a e ^(M^") we define the Weyl quantization 
of the symbol o, the operator 

Op(a) = / a(T],3;)ei(^-^+^-^) (157) 

whose action on any X]/ e L^(R") is 

{Op{a)xi/){x) = [ Ka{x,y)xi/{y)dy (158) 

with a kernel 

Observe that the kernel Ka defined in (159) is obtained from <7 by a partial Fourier 
transform followed by an invertible and measure-preserving change of variables. These 
operations make sense when o is an arbitrary temperated distribution and define Ka as 



a temperated distribution. Therefore, it is possible to define the quantization of a generic 
temperated distribution. Later we will use the quantization of the Dirac 5 distribution. 

Now we want the inverse mapping of the Weyl quantization, namely we want to 
obtain the symbol starting from the operator. From (159) it is easy to check that the 
dequantization of the operator S whose kernel is K{x,y) is given by the symbol 

G{q,p)^Or\S){q,p) = {2nhrW{K){q,p), (160) 

where 

»'(^)fe/')=X/(«+i.*-f)e-'"'»(^ (160 

is known as the Wigner transform of F e L?-(E?"). This is the fundamental connection 
between Weyl quantization and Wigner transform. By a straightforward direct computa- 
tion, that we leave as an exercise, one has the following remarkable relation: 

Tr(0p(a)0p(i:)):= / Ka{x,y)K,{y,x) dxdy ^ [ c>{q,p)z{q,p) (162) 

for any a, T G ^(M^"). This property will be crucial for our tomographic framework. 

The statistical states of a quantum system are described by positive (self-adjoint) 
operators with unit trace on the Hilbert space L^(]R"), namely p >0, p* = p and Trp = 1, 
that are called density operators. It is easy to check that a density operator p is an integral 
operator on L^(M"), and if 

p ^Y^^kVkiVk,-) (163) 
k 

is the spectral resolution of p, with { V^} an orthonormal eigenbasis and 1 > Ay^ > Xk+i > 
the corresponding eigenvalues, the trace condition is equivalent to Y,k^k = 1> ^iid the 
kernel of p is (with a little abuse of notation) 



Pix.y) = Y,^k¥kix)\i/kiy)- (164) 

k 

The expectation value in the state p of an observable associated to a bounded operator S 
is given by 

(S)p = Tr(5p) = Y.Xk{Vk,SWk)- (165) 
k 

Therefore, the density matrix p is given by a convex combination of wave functions Xj/^ 

with weights Xk- Physically, it represents a classical statistical mixture of pure quantum 
states extracted with probabilities X^. By making use of Eqs. (160) and (162), the 
quantum expectation value (165) of 5 = Op (a) reads 

(Op(c7))p=Tr(Op((7)p)= / G{q,p)Wip){q,p)dqdp= {a)wip). (166) 



This equation expresses the expectation value of a quantum observable in the quantum 
state p as a classical expectation value of its symbol in the "classical" state W{p). In 



this respect, W{p) is called Wigner quasi-distribution function of p and is a substitute 
for the nonexistent joint probability density of momentum and position in the quantum 
state p. Indeed, since the uncertainty principle imposes a limit on the precision with 
which momentum and position can be determined at the same time in the state p, it does 
not make sense to speak of a joint probability distribution for these observables. In fact, 
W{p) is not a genuine probability density, because it may assume negative values (e.g. 
if p = i^(v/,-), with V/ G odd, we have that W(p)(0,0) = -\\\if\\l/{nhY < 0). 

It is, however, always real and in some sense it tries very hard to be a joint density for 
momentum and position. The above discussion provides supporting evidence for this 
heuristic assertion. Moreover, it is easy to check that W{p) has the right marginals: 

I W{p)iq,p)dp = p{q,q) = Y.h\¥k{q)\^- (167) 

Now we can apply all the results obtained in Sees. 3 and 5 to reconstruct the state of 
a quantum system using a tomographic approach. In the quantum framework the object 
that we want to reconstruct is the state p of the system. Therefore, let us compute the 
expectation of 

5(X-fl-q-V-p) ■.= Op{5{X-fl-q-V-p)) (168) 
in the state p. By (166) and (151) one has 

{5{X-fl-q-V-p))p = {d{X-fl-q-y-p))^^p)=W{pf\x,fl,v). (169) 

Therefore, by measuring the expectation in the state p of the family of operator-valued 
distributions (168), with X eR, (n^v) e M^", we in fact are gathering all the tomograms 
of the Wigner function W^(p). Therefore, by applying the inversion formula (151) we can 
recover W{p), that is, up to a normalization factor, the symbol of the density operator p 
[see Eq. (160)]. Thus, by quantizing it we can reconstruct the quantum state p. Explicitly, 
by using Eqs. (160), (159) and (151), we get 

pix,y) = _^^W(p)(^,77)ei''(--^)/^di7 

= / dr]e^<^~yy^[ ^^^w(p)^'(X,Ai,v)eH^-^-(-+^)/2---^) 

= j W{pf\x,^,v)ei^-^<^^yy')5{v-ix-y)/h)^^ 

= f W(p)^^Z,M,(x-3;)A)ei(^-^-(-+^)/2)^. (170) 

We recall that the tomogram W(p)^^^^ is the marginal probabiUty density of the quan- 
tum observable Xi^^y^{q,p) = fi-q + v-p. When = (jU, v) e S^""^ the observable 



is nothing but the component of a generalized position, called quadrature, in a ro- 
tated frame of the phase space. From the experimental point of view, quadratures of 
photons can be measured in a typical quantum optics setup by a homodyne detection 
scheme. A different scheme was proposed for the longitudinal motion of neutrons [3]. 
Recently, in a spectacular experiment [13] quadratures were measured for reconstructing 
the transversal motional states of helium atoms, by using a free evolution of the wave 
packet followed by a position sensitive measurement. 

The determination of the quantum state represents a highly nontrivial problem, whose 
history can be traced back to the early days of quantum mechanics, namely to the 
Pauli problem of reconstructing a wave function from its marginals in position and 
momentum. The experimental validation had to wait until quantum optics opened a new 
era. In this section we have followed the idea of Vogel and Risken [24]. Since then 
many improvements and new techniques have been proposed [15, 11, 19]. An up-to- 
date overview can be found in Ref. [21]. 



Notation 





Radon transform and its dual, see Eq. (12) and Eq. (13) 


A 

h 


Tomogram of / along the direction t, , see Eq. (23) 




transform, see Eq. (72) 


y{M) 


Schwartz space of rapidly decreasing functions on M, see Eq. (16) 


5{x) 


Dirac delta distribution at 


fj 


Fourier transform and its inverse, see Eqs. (37)-(38) 


W{p) 


Wigner transform of the state p, see Eq. (161) 


Op(C7) 


Weyl quantization of the symbol <J, see Eq. (157) 


pn 


Space of all the {n - J) -dimensional affine subspaces of M". (P" = P") 


K 


= M"\{0}, M+ = [0,+oo), N*=M\{0} 




= / A^^^ e""^ dA , Euler gamma function 



^0 
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